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We use time- and angle-resolved photoemission to measure quasiparticle relaxation dynamics 
across a laser-induced superconducting phase transition in BES^CaCi^Os+s- Whereas low-fluence 
measurements reveal picosecond dynamics, sharp femtosecond dynamics emerge at higher fluence. 
Analyses of data as a function of energy, momentum, and doping indicate that the closure of the 
near-nodal gap and disruption of macroscopic coherence are primary mechanisms driving this onset. 
The results demonstrate the important influence of transient electronic structure on relaxation 
dynamics, which is relevant for developing an understanding of nonequilibrium phase transitions. 


In the study of complex many-body interactions, a fun¬ 
damental topic is the behavior of materials at the bound¬ 
aries of ordered phases, which can be accessed, for ex¬ 
ample, by monitoring material properties as a function 
of temperature, magnetic field, or chemical doping. Re¬ 
cently, there has been great interest in studying phase 
transitions not just in steady state, but also out of equi¬ 
librium [1, 2]. It has been shown, for example, that ul¬ 
trafast optical pulses can probe and in cases even ma¬ 
nipulate electronic order in materials ranging from high- 
temperature superconductors [3-10], to topological insu¬ 
lators [11, 12], to charge-density-wave materials [13-15]. 
Although increasingly powerful theoretical techniques are 
being developed [16-20], our understanding of these phe¬ 
nomena remains lacking. 

In high-temperature superconductors, particularly in 
the cuprates, an important related question is how 
quasiparticle interactions are influenced by electronic 
structure when the system transitions across the phase 
boundary between the superconducting and pseudogap 
regimes. An intriguing recent approach for insight into 
this physics has been the examination of quasiparticle 
relaxation using pump-probe spectroscopy at low tem¬ 
perature but high pump fluence, where the supercon¬ 
ducting condensate is fully vaporized into quasiparti¬ 
cles and must dynamically re-emerge. This regime has 
been explored using time-resolved reflectivity and trans¬ 
missivity [3, 21-23], and studies have reported distinct 
femtosecond and picosecond relaxation time scales that 
were often assigned to the condensate and pseudogap 
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phases. However, such all-optical probes lack momen¬ 
tum resolution, and have no direct access to the dynam¬ 
ics of the electronic band structure. Conversely, time- 
and angle-resolved photoemission spectroscopy (time- 
resolved ARPES) [24-27] is ideally suited for probing 
both momentum-dependent quasiparticle [8, 28-35] and 
gap [8, 36, 37] dynamics in cuprates, yet low-temperature 
studies using fluences high enough to destroy supercon¬ 
ductivity [30-32, 34-37] have to date provided no com¬ 
parisons of the gap and quasiparticle population. 

In the present study, we use time-resolved ARPES to 
probe the critical-fluence regime of the high-temperature 
superconductor Bi 2 Sr 2 CaCu 2 08+<5 (Bi2212) in detail, 
providing a characterization of quasiparticle relaxation 
rates as a simultaneous function of crystal momentum 
and energy in cuprates. Fluences inducing a complete 
closure of the near-nodal gap result in a distinctive 
two-component relaxation signature, with well-defined 
femtosecond and picosecond recovery scales that occur 
throughout a large range of k- space. Correlations be¬ 
tween the onset of these two components and gap dy¬ 
namics establish a strong connection between picosec¬ 
ond relaxation dynamics and superconductivity. How¬ 
ever, the detailed momentum and doping dependence of 
the faster, femtosecond component contradicts an assign¬ 
ment to pseudogap-phase dynamics. Instead, the dy¬ 
namics can be explained by release of kinetic restric¬ 
tions within a “normal-state” phase obtained by a full 
quenching of the condensate. The results demonstrate 
the important influence of transient electronic structure 
on relaxation dynamics, which is relevant for developing 
an understanding of nonequilibrium phase transitions. 

The experimental apparatus is as described in Ref. [27] , 
and uses 836-nm (hv = 1.48 eV) pump pulses and 209-nm 
(hu = 5.93 eV) probe pulses. Energy, momentum, and 
time resolutions are 23 meV, 0.003 A -1 , and 300 fs. We 
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FIG. 1. Evolution of the quasiparticle population and band 
gap in superconducting Bi2212 (T c = 91 K) following an ul¬ 
trafast infrared pump pulse, (a)—(b) Equilibrium ARPES 
dispersions (t — —1.1 ps) for cuts at 4> = 45° and 0 = 30° (see 
panel (a) schematic), (c)— (d) Nodal (0 = 45°) quasiparticle 
population SI(t), from integrating ARPES intensity between 
the white double arrows in (a). Main panels show T — 20 
K. Inset shows data at 23 /iJ/cm 2 for T = 100 K (T > T c ). 
(e)— (f) Normalized superconducting gap A (t) from the cut 
at (j) — 30°, extracted by fitting symmetrized EDCs at fcp to 
a broadened BCS line shape [36]. 


have corrected for detector nonlinearity [27], as well as 
for a pump-induced time-dependent but uniform energy 
shift in the electronic spectrum (< 4 meV) [36, 38]. 

Figure 1 shows a fluence-dependent analysis of a nearly 
optimally doped sample (T c = 91 K), where two cuts 
through k -space are depicted, at 0 = 45° and 0 = 30° 
(0 is defined from the Y point relative to Y-M as in 
the panel (a) inset). The cut at 0 = 45° intersects a 
d-wave gap node, where quasiparticle dynamics can be 
measured independently of gap dynamics, and can be 
characterized with a high signal-to-noise ratio (Figs. 1(c) 
and 1(d)) by integrating spectral intensity change across 
a large window in energy and momentum (white dou¬ 
ble arrows in Fig. 1(a)). Concomitant gap dynamics 
(Figs. 1(e) and 1(f)) are extracted using energy distri¬ 
bution curve (EDC) symmetrization [8, 36, 39] from the 
cut at (j) = 30°, which has an equilibrium gap of 14 meV 
(Fig. 1(b)). Consistent with previous measurements [8], 
at lower fluence the gap magnitude is suppressed but al¬ 
ways finite, and the population of nonequilibrium quasi¬ 
particles decays with picosecond dynamics. As shown in 
Figs. 1(d) and 1(f), however, beyond a critical fluence F c 
(defined as the fluence necessary to close the near-nodal 
gap) [36, 37], a sharp femtosecond relaxation signature 
emerges in the quasiparticle population response. The 
signature is also coupled to the near-nodal gap tempo¬ 
rally, appearing most prominently when the gap is closed 


FIG. 2. (a) Amplitudes SIosiow and Slofast , from fits to the 

data in Figs. 1(c) and 1(d) as described in the text [43]. (b) 
Minimal gap from Figs. 1(e) and 1(f) versus fluence. (c)— 
(d) Time constants r s i ow and r/ as t, from fits to the data in 
Figs. 1(c) and 1(d) as described in the text [43]. 


(t < 1.5 ps, see shaded gray regions in Figs. 1(d) and 
1(f)). We note that, although two-component quasipar¬ 
ticle recovery dynamics have been previously observed in 
cuprate studies using time-resolved ARPES [29-31] and 
time-resolved reflectivity and transmissivity [3, 4, 21- 
23, 40-42], this connection between two-component be¬ 
havior and the superconducting gap has not been re¬ 
ported. 

Figure 2 shows a more detailed analysis of quasiparticle 
decay, where fast and slow component amplitudes and 
time constants are extracted from Fig. 1 by fitting the 
integrated ARPES intensity with a bi-exponential decay 


SI{t) = ©(£) 


Cfaste - tlT + C^e-^ 


( 1 ) 


after convolution with a Gaussian of 300-fs duration to 
incorporate time resolution. In the above, t = Delay —to, 
and t f as t and r s i ow are the time constants of the respec¬ 
tive fast and slow components of the relaxation, while 
Cfast and C s i ow are their bare amplitudes. After convo¬ 
lution with the resolution function, we extract the effec¬ 
tive amplitudes Slofast and SIo s i ow , which directly relate 
to the data [43]. For the 4 fi J/cm 2 measurement we con¬ 
strained Cfast = 0 because a fast component was not 
apparent. We also measured the above-T c response, set¬ 
ting C s i ow = 0 and constraining t < 0.6 ps because a 
slow component was not apparent. 

Overall, the picosecond quasiparticle recovery compo¬ 
nent exhibits many trends that establish a close connec¬ 
tion to superconductivity. For example, the picosecond 
component only appears in the low-temperature data. As 
shown in Fig. 2(c), when present, r s i ow decreases with 
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FIG. 3. Energy and momentum dependence of quasiparti¬ 
cle recombination dynamics in optimally doped Bi2212 (25 
/xJ/cm 2 pump fluence, T — 20 K, T c = 91 K). (a)-(b) 
Energy-resolved ARPES intensity change 51 (t), momentum- 
integrated between \k — &f| < 0.08 7r/a for nodal, off-nodal, 
and far-off-nodal cuts, at two representative energies above 
Ef . (c)—(d) Energy-resolved quasiparticle relaxation time 

constants and amplitudes. The red, blue, and green vertical 
dashed lines mark the gap edge at 0 = 45°, (f> = 37°, and 
(p = 26°. Intensities in (d) are obtained by scaling momen¬ 
tum cuts so that f_ Q e 3eV I(k,u) du averages to 1 for k such 
that — 0.0757r/a < k — kF < —0.037T /a when t < 0. 


increasing fluence [8], as expected of quasiparticles reen¬ 
tering a superconducting condensate following rules of 
second-order kinetics [5, 8, 44]. The amplitude 5Iq s i ow 
also increases rapidly with fluence in the low-fluence 
limit, yet appears to saturate above F c (Fig. 2(a)), which 
is consistent with a density of quasiparticles decaying 
into the superconducting state in a manner fundamen¬ 
tally limited by the equilibrium superfluid density. By 
contrast, the femtosecond recovery component exhibits 
very different behavior. At low temperature, 5Iof ast is 
suppressed or absent at low fluence and only becomes 
substantial for F > F c . At high temperature (above 
T c ), femtosecond dynamics dominate. At all tempera¬ 
tures, t f as t exhibits no fluence dependence or may even 
increase with fluence (Fig. 2(d)). 

Figure 3 shows quasiparticle relaxation dynamics, re¬ 
solved in both energy and momentum, for momentum 
cuts at cf) = 26° (outside the pseudogap-state Fermi arc), 
(j) = 37° (inside the pseudogap-state Fermi arc), and 
< j) = 45° (the node), responding to a pump fluence of 


25 fi J/cm 2 (above F c ). Such an analysis takes advan¬ 
tage of the full power of time-resolved ARPES, and two- 
component relaxation dynamics notably appear through¬ 
out a large portion of the Brillouin zone. As clear from 
Figs. 3(a) and 3(c), quasiparticle relaxation after the first 
picosecond occurs significantly faster away from the node 
than at the node, in agreement with previous energy- 
integrated measurements at lower fluence [8] . These mo¬ 
mentum variations in r s i ow become larger with decreasing 
energy (Fig. 3(b)), and establish further connections to 
superconductivity, as low-energy quasiparticles necessar¬ 
ily interact more strongly with the many-body ground 
state than do their higher-energy counterparts [45]. 

Though it is tempting, based on Figs. 1 and 2, to as¬ 
sociate the fast component with a competing pseudogap 
order, Fig. 3 shows that the momentum- and energy- 
dependent trends in 5Iof ast are directly at odds with 
the fact that the pseudogap becomes most prominent to¬ 
ward the Brillouin zone face. At optimal doping, the 
equilibrium pseudogap is absent or undetectable near 
4> = 45°, and becomes increasingly pronounced toward 
4> = 0°, starting at « 30° [47, 48]. As shown in 
Fig. 3(a), by contrast, at low energy the fast component 
completely disappears in the off-nodal data. In part, this 
reflects the fact that the low-energy off-nodal quasiparti¬ 
cle relaxation signature is overshadowed by the dynam¬ 
ics of the nonequilibrium superconducting gap. How¬ 
ever, even above the gap edge, Tf ast exhibits no dis¬ 
cernible momentum dependence within experimental er¬ 
ror (Figs. 3(b) and 3(c)), and energy-dependent ampli¬ 
tudes 5Iof as t do not appear to increase near the Brillouin 
zone face (Fig. 3(d)). Rather, they may even tend to de¬ 
crease, although ARPES matrix elements make it some¬ 
what difficult to make these comparisons quantitative. 

Stronger evidence differentiating the fast component 
from the physics of the pseudogap occurs in the doping 
dependence of quasiparticle relaxation. Figure 4 shows 
nodal quasiparticle relaxation dynamics at comparable 
excitation densities (« 24 /iJ/cm 2 ) for multiple dop¬ 
ings of Bi2212, corresponding to critical temperatures 
T c = 78 K (underdoped, UD78K), T c — 91 K (nearly op¬ 
timally doped, OP91K), T c = 78 K (overdoped, OD78K), 
and T c = 59 K (very overdoped, OD59K). For the 
first three dopings, both equilibrium ARPES [46, 47] 
and time-resolved ARPES [32] measurements have previ¬ 
ously identified distinct transition temperatures T c and 
T* for the onset of superconductivity and the pseudo¬ 
gap, with T* occurring at 200 K for the UD78K sam¬ 
ple, 150 K for the OP91K sample, and 97 K for the 
OD78K sample. The OD59K sample has no equilibrium 
pseudogap [46, 47]. In spite of these differences, two- 
component recovery dynamics prominently appear in the 
time-resolved nodal quasiparticle relaxation dynamics for 
all four dopings (Figs. 4(a)-4(d)), and quasiparticle re¬ 
laxation above T c in the OD59K sample resembles the 
low-temperature femtosecond component for the OD59K 
sample (Fig. 4(d)) similarly to the way that quasiparticle 
relaxation above T c , but below T*, in the OP91K sam- 


























4 



FIG. 4. Doping dependence of quasiparticle decay for equilibrium temperature T — 20 K (unless otherwise noted) and fluence 
F — 23 /iJ/cm 2 . (a)—(d) Nodal quasiparticle response curves, extracted as in Fig. 1. The fits are bi-exponential decay functions 
(Eq. 1). Panel (d) inset shows the OD59K nodal quasiparticle response at T = 70 K. (e)—(g) Cartoon illustrations of the 
impact on quasiparticle scattering of opening a d -wave gap. (h)-(i) Time constants from the fits in (a)-(d) [43]. (j) Amplitude 
component ratios SIo slow/(bio slow + 8Io fast) from the fits in (a)-(d). Inset is the equilibrium Bi2212 phase diagram, showing 
the relationship between doping, the pseudogap (PG), and superconductivity (SC) [46, 47]. 


pie resembles the OP91K low-temperature femtosecond 
relaxation component (Fig. 1(d)) [43]. 

Taking the momentum and doping dependence of 
quasiparticle relaxation rates into consideration jointly, 
we argue that the primary explanation for the onset of 
two-component dynamics is that quasiparticle relaxation 
dynamics are directly affected by changes that occur in 
the quasiparticle spectrum as the system relaxes back to¬ 
ward equilibrium. More specifically, the opening of a gap 
dramatically reduces the amount of available phase space 
near the Fermi level. A quasiparticle in metallic Bi2212 
can relax by many channels (Fig. 4(e)), with a final-state 
phase-space density roughly proportional to the product 
of the Fermi surface area and the quasiparticle energy 
(assuming a uniform Fermi velocity). In the presence of 
a d-wave gap, however, scattering rates for lower-energy 
states are sharply curtailed (Fig. 4(f)). For a quasipar¬ 
ticle residing along the nodal direction at 15 meV in 
optimally doped Bi2212 (A 0 ~ 35 meV), the gap im¬ 
poses an almost 80% phase-space reduction for relaxation 
via scattering (particle-conserving) interactions [49] . Be¬ 
yond this, the gap opens in the midst of the quasiparticle 
relaxation process, which means that quasiparticles are 
steadily lifted from lower energies back to higher energies 
on a picosecond timescale (Fig. 4(g)), and the overall ef¬ 
fect counterbalances quasiparticle relaxation at fixed en¬ 
ergy. This explains, for example, why the energy of the 
maximum r s i ow increases with decreasing (j) in Fig. 3(c). 
While the effect occurs most dramatically at the Bril- 
louin zone face, nodal and near-nodal states might also 


be impacted because of the resultant conversion of boson- 
absorption scattering channels into boson-emission chan¬ 
nels. We note that momentum conservation requirements 
are ignored in these first two arguments, which may 
be reasonable given the large inhomogeneities that are 
known to exist in Bi2212 on nanometer length scales [50]. 
Finally, there is a fundamental change—due to coher¬ 
ence factors—in electron-boson coupling matrix elements 
in the presence as opposed to the absence of supercon¬ 
ductivity, which may affect quasiparticle recombination 
rates in important ways [45, 51]. For example, if the 
dominant channel for quasiparticle recombination is be¬ 
tween k -space regions where the gap parameter has 
like signs, and the dominant scattering interaction is odd 
under time reversal, then recombination will be reduced 
in the superconducting state as compared to the normal 
state. 

Curiously, Figs. 4(h) and 4(i) reveal that time con¬ 
stants for both the fast and slow components decrease 
sharply with decreasing hole concentration, and the am¬ 
plitude ratio 5Io s i ow /SIo total is markedly different for the 
OD59K sample than it is for the other three dopings 
(Fig. 4(j)). Further studies with better time resolution 
are needed to investigate these effects more fully. 

In summary, we have shown that kinetically released 
quasiparticle relaxation dynamics and the opening of the 
superconducting gap play a large role in generating two- 
component relaxation dynamics, independently of the 
role that the pseudogap might play. This result will be 
important for any optically-induced phase transition in- 
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volving a gap, but is particularly relevant for cuprates 
given the findings of many all-optical studies, which also 
report a femtosecond relaxation component, yet have as¬ 
sociated it in many cases as a direct manifestation of the 
pseudogap [3, 4, 21-23, 40-42]. We note that all-optical 
studies are typically sensitive to energies far higher than 
those probed in the current work. Hence, the two tech¬ 
niques may be accessing different dynamics if, for ex¬ 
ample, the influence of the pseudogap extends to higher 
energies than does the superconducting gap, or if the 
low-energy signatures of the pseudogap are more local¬ 
ized in momentum than higher-energy signatures. Time- 
resolved ARPES studies at higher fluences or closer to 


the first Brillouin zone boundary may be helpful in re¬ 
solving this question unambiguously. 
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